Diquark condensation in dense SU(2) matter 
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We report on a lattice study of two-color QCD with adjoint staggered fermions at high density. We find that the 
model has no early onset and we report on results for diquark condensation, from simulations with and without 
a diquark source term. 



1. INTRODUCTION 

Numerous model calculations suggest that the 
ground state of QCD at high baryonic density 
is characterized by diquark condensation which 
spontaneously breaks gauge and/or baryon num- 
ber symmetry (see jl],U for recent reviews). In 
|H^] we presented evidence that two-color QCD 
with one flavor of adjoint staggered fermions has 
interesting QCD-like features, in particular the 
absence of baryonic Goldstone modes in the spec- 
trum. Here we extend the study of this model to 
include possible diquark condensates. We refer 
to H for a detailed description of the model and 
its properties, while we summarize here only the 
main features: N flavors of staggered fermion, 
in the adjoint representation of the SU(2) gauge 
group, have a XJ(2N) chiral-flavor symmetry (at 
zero mass m and zero chemical potential /i) . This 
is explicitly (spontaneously) broken by the mass 
(the chiral condensate) to Sp(2A). SSB gives rise 
to 2A 2 - 1 Goldstone modes, of which N(N - 1) 
are baryonically charged. When N > 1 and 
fi > the vacuum starts rotating from the 

chiral condensate direction to that of a diquark 
condensate ||, which in the case of N = 2 is 
(modulo a baryonic rotation): 




where e p<3 , p, q = 0,1 is the completely anti- 
symmetric tensor in flavor space. If N = 1 no 



baryonic Goldstone mode is expected, and no lo- 
cal, gauge invariant, Lorentz scalar diquark con- 
densate is possible. Therefore no early onset 
(/i ~ m^) transition is expected. One interesting 
possibility is a gauge non singlet — hence color 
superconducting — diquark condensate: 

<?<& = \ [X tr {x)t l x{x) + x(x)t l X tr (x)] , 

where t 1 , i = 1, 3 are color generators. Note how- 
ever that the fermionic determinant of the N = 1 
model (although real) does not have a definite 
sign. The indication in [||JJ] was that the effect 
of the sign is indeed important. If we ignore the 
sign, the model is well described by Chiral Pertur- 
bation Theory (xPT) ||, including the presence 
of baryonic Goldstone modes. The latter are sup- 
pressed only by the inclusion of the sign. Here we 
want to complement the previous analysis with 
the study of diquark condensation ||. 

2. ALGORITHMS 

We studied the N — 1 model by means of both 
HMC § and TSMB § algorithms. Because of 
lack of ergodicity, with HMC we could explore 
only the positive determinant phase. TSMB is 
defined by a set of parameters (degrees of poly- 
nomia n^, the number of heatbath and overre- 
laxation iterations for the bosonic fields It and 
I Q , the number of metropolis iterations for the 
gauge fields Im and the number of noisy correc- 
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Nf^ q is the autocorrelation length (expressed in 
number of matrix multiplications) for the plaque- 
tte, evaluated from a run of N sw sweeps, n 2 = 90 
at /i = 0, otherwise n 2 ~ 10ni. 
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N g is the effective number of independent config- 
urations, is the fraction with a negative de- 
terminant and (r) the average reweighting factor 
including the sign. 



tions I c ). One has to tune such parameters in or- 
der to minimize the autocorrelation. The latter 
must be expressed in term of matrix multiplica- 
tions. The number of matrix multiplications per 
sweep is roughly given by: 

7 

Nmuit ^ -^ni{I h + I + I M ) + {n 2 + n 3 )I c . 

In general when the chemical potential increases 
the lowest eigenvalues (e) approach zero and the 
simulations become difficult. In Table [l] we show 
how the condition number (A/e, A being the 
largest eigenvalue) as well as the plaquette au- 
tocorrelation change by increasing pL. The au- 
tocorrelation of fermionic observables is typically 
about 10 times shorter. From Table || we see that 
until the point /i = 0.38 the average sign is still 
not too small, while points above pL = 0.4 are 



already exceedingly difficult. 

3. DIQUARK CONDENSATES 

Our model has very interesting features, which 
distinguish it from other SU(2) models. It also 
has a sign problem. It is natural and interesting 
to study the effect of the inclusion of the sign. 
We distinguish two models: the one where only 
positive (determinant) configurations are sam- 
pled and the one where both positive and neg- 
ative configurations are included with their sign. 
The positive sector turns out to be in very good 
agreement with %PT predictions; i.e., the positive 
model behaves effectively like an N — 2 model. 

In order to study possible diquark condensa- 
tions we introduced a source j for qq$. This is 
done in a partially quenched way: by introducing 
the source only in the measurements. xPT pre- 
dicts also the dependence on j of the chiral and 
the diquark condensate || . In particular the sum 

(^)o = W) 2 + (qqsf (1) 

must, to lowest order in x?T, be a constant inde- 
pendent of both on n and j (for NLO see Q ) . In 
fact, as one can see on fig. [I], this relation is very 
well satisfied in the positive sector. The validity 
of the relation (|l|) also allows an extrapolation of 
{qqs) to zero j, which would otherwise be quite 
difficult. 

Since the condensate qq sc is not gauge invariant 
we studied its susceptibility Xsc 

Xsc = \(x{x)t l x tr {x)x tr {xYx{x)). 

This has a finite value even without a source term. 
In the low density regime \sc shows little varia- 
tion with fi. Once we reach the X-PT transition 
point, Xsc drops rapidly in the positive sector. 
We interprete this as an effect of Pauli blocking 
or phase space suppression: as the ground state 
is filled up by fermions there is less phase space 
left to accommodate the fermion loops that con- 
tribute to Xsc- When the sign is properly taken 
into account this effect disappears and Xsc re- 
mains stable through the transition (fig. ||). 
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Figure 1. Dotted (increasing) lines are qq3, 
dashed (decreasing) lines are (^ip) (at different 
values of \i between and 0.5). Solid lines are 
the sum in (|l|). The error bars do not take into 
account the autocorrelation. 



4. CONCLUSIONS 

The interesting lesson one can learn from this 
model is that the early onset transition, typical of 
all other SU(2) models, is indeed delayed as ex- 
pected from the analysis of possible diquark con- 
densates. Numerically this is obtained by a deli- 
cate balance between configurations with positive 
and negative determinant. This makes simula- 
tions prohibitively expensive in the region where 
one would expect a true onset transition at high 
density. 
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